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Level densities for ^^~^*Mo have been extracted using the ('^He,Q7) and ('^He,'^He'7) reactions. 
From the level densities thermodynamical quantities such as temperature and heat capacity can be 
deduced. Data have been analyzed by utilizing both the microcanonical and the canonical ensemble. 
Structures in the microcanonical temperature are consistent with the breaking of nucleon Cooper 
pairs. The S-shape of the heat capacity curves found within the canonical ensemble is interpreted 
as consistent with a pairing phase transition with a critical temperature for the quenching of pairing 
correlations at Tc ~ 0.7 - 1.0 MeV. 
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I. INTRODUCTION 

Level density is a characteristic property of many-body 
quantum mechanical systems. Its precise knowledge is 
often a key ingredient in the calculation of different pro- 
cesses, such as compound nuclear decay rates, yields of 
evaporation residues to populate exotic nuclei, or ther- 
monuclear rates in astrophysical processes. 

Measurements of experimental nuclear level density are 
an important prerequisite for thermodynamical studies of 
atomic nuclei. Level density is directly connected to the 
multiplicity of states, i.e., the number of physical realiza- 
tions of the system at a certain excitation energy. The 
entropy is a fundamental quantity and a measure of the 
disorder of the many-body system. Within the micro- 
canonical ensemble it is defined as the natural logarithm 
of the multiplicity of states. When the entropy is known, 
thermodynamic quantities such as temperature and heat 
capacity can be extracted. These quantities depend on 
the statistical properties of the nuclear many-body sys- 
tem and may reveal phase transitions. 

Pairing correlations are one of the fundamental proper- 
ties of nuclei, and have been successfully described by the 
Bardeen-Cooper-Schrieffer (BCS) theory of superconduc- 
tivity 01 . By using the BCS theory the thermodynami- 
cal properties of nuclear pairing were investigated in the 
study of warm nuclei Q, M Qi IB| - In the case of a finite 
Fermi system such as the nucleus, statistical fluctuations 
beyond the mean field become important. The fluctua- 
tions smooth out the sharp phase transition, and then the 
pairing correlations do not vanish suddenly but decrease 
with increasing temperature. The quenching of pair- 
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ing correlations has been obtained in recent theoretical 
approaches: the shell model Monte Carlo (SMMC) cal- 
culations 0, S 13 ) the finite-temperature Hartree-Fock- 
Bogoliubov theory , and the relativistic mean field the- 
ory Experimental data on the quenching of pair 
correlations are important as a test for nuclear theories. 
A long-standing problem in experimental nuclear physics 
has been to observe the transition from strongly paired 
states at zero temperature to unpaired states at higher 
temperatures. A signature of the pairing transition at fi- 
nite temperature might be a local increase in the heat 
capacity as a function of temperature [Tlj . Recently 
[12. fine structures in the level densities in the 1- 
7 MeV region were reported, which are probably due to 
the breaking of individual nucleon pairs and a gradual 
decrease of pair correlations. 

The group at the Oslo Cyclotron Laboratory (OCL) 
has developed a method to extract simultaneously the 
level density and the radiative strength function from 
primary 7 spectra IT^ . The method is a further de- 
velopment of the sequential extraction method [TH IT^ . 
The Oslo method has been tested in the rare-earth 
mass re gion which led to many interesting applications 
[T^ Ir^ list figf . In order to make quantitative judg- 
ments of the applicability of the method, the Oslo Cy- 
clotron group has extracted the level density and radia- 
tive strength function (RSF) of the very light ^^'^^Si nu- 
clei, where these quantities are known. Excellent overall 
agreement was found [23|. Subsequently, another exten- 
sion has been made to the intermediate nuclei ^^'^'^Fe and 
96,97]yjQ^ and it has been shown that the method can be 
applied in this intermediate mass region where the level 
density is still relatively low [2lll2^ . All of these success- 
ful applications have motivated us to employ the Oslo 
method to study medium-heavy nuclei in the vicinity of 
closed shells. 
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The naturally occurring isotopes of molybdenum span 
one of the larger isotopic ranges and are well suited as tar- 
gets for the study of nuclear properties, such as the effect 
of changing from spherical to deformed shapes. When 
approaching closed shells, the nuclear structure changes 
significantly, and one expects this to influence the level 
densities and radiative strength functions. 

The even-even ^^Mo has a filled TV = 50 neutron shell 
[23| . It is essentially a spherical nucleus and vibrations 
are primarily governed by the proton core. As the mass 
increases from ^''Mo to ^'^°Mo, neutrons fill the 2^5/2 and 
I57/2 sub-shells. Moving away from the iV = 50 shell clo- 
sure, pairing and quadrupole interactions cause a more 
collective behavior in the heavier Mo isotopes. The char- 
acter of the isotopes changes rapidly from that of the 
essentially spherical ^^Mo to nuclei making a transition 
from collective vibrators to the deformed rotors of the un- 
stable i°'*Mo and ^°^Mo isotopes The transitional 
nature of molybdenum isotopes away from TV = 50 has 
been the focus of several efforts as described in |2^ and 
references therein. 

Around closed shells, effects from the increasing single 
particle energy spacings can be expected. These will also 
influence the entropy difference between odd-mass and 
even-even nuclei. Therefore, a statistical description of 
the transition from closed shells to deformed nuclei is of 
great interest. 

In this work, a unique and consistent investigation of 
the six ^•^~^^Mo isotopes is performed in order to deter- 
mine experimentally the level density from the ground 
state up to the neutron binding energy. The Oslo method 
also determines the RSFs of the molybdenum isotopes 
studied; these are presented in an earlier paper . 



II. EXPERIMENTAL METHODS 

The experiments were carried out at the Oslo Cy- 
clotron Laboratory by bombarding 9^,96, 97, gs^yj^ targets 
with ^He ions. In the present work, results from eight 
different reactions on four different targets are discussed. 
These are the following six reactions which are the sub- 
ject of the present investigation: 

1. 98Mo(3He,a7)97Mo (45 MeV) 

2. 98Mo(3He,3He'7)9SMo (45 MeV) 

3. 96Mo(3He,a7)95j4o (30 MeV) 

4. 96Mo(3He,3He'7)96Mo (30 MeV) 

5. 94Mo(3He,a7)93]y[Q ^eV) 

6. 94Mo(3He,3He'7)94Mo (30 MeV) 
together with the reactions 

7. 97Mo(3He,a7)9^Mo (45 MeV) 

8. 9^Mo(3He,3He'7)9^Mo (45 MeV) 

which have been reported earlier [21I |2^ . The self- 
supporting targets with thicknesses of ^ 2 mg/cm^ are 
enriched to ~ 95%. The experiments were run with beam 
currents of ~ 2 nA for 1-2 weeks. The particle-7 coinci- 
dences were measured with the CACTUS multi-detector 
array. The charged ejectiles were detected by eight par- 



ticle telescopes placed at an angle of 45° relative to the 
beam direction. An array of 28 Nal 7-ray detectors with 
a total efficiency of ^15% surrounded the target and par- 
ticle detectors. 

For each initial excitation energy, the 7-ray spectra 
are recorded as a function of the initial excitation energy 
of the residual nucleus. This is accomplished by utiliz- 
ing the known reaction Q-values and kinematics. Using 
the particle-7 coincidence technique, each 7 ray can be 
assigned to a cascade depopulating a certain initial ex- 
citation energy in the residual nucleus. The data are 
therefore sorted into total 7-ray spectra originating from 
different initial excitation-energy bins. Each spectrum 
is then unfolded with the Nal response function using a 
Compton-subtraction method which preserves the fluc- 
tuations in the original spectra and does not introduce 
further, spurious fluctuations [27| . From the unfolded 
spectra, a primary-7 matrix P{E^ Ey ) is constructed us- 
ing the subtraction method of Ref . |23| . 

The basic assumption of this method is that the 7-ray 
energy distribution from any excitation energy bin is in- 
dependent of whether states in this bin are populated 
directly via the (■^He,a) or (^He,'^He') reactions or indi- 
rectly via 7 decay from higher excited levels following 
the initial nuclear reaction. This assumption is trivially 
fulfilled if one populates the same levels with the same 
weights within any excitation energy bin, since the decay 
branchings are properties of the levels and do not de- 
pend on the population mechanisms. The assumptions 
behind this method have been tested extensively by the 
Oslo group and have been shown to work reasonably well 

The (^He,'^He'7) and (■^He,a7) reactions exhibit very 
different reaction mechanisms. This is demonstrated in 
Fig. n where the particle spectra in coincidence with 7 
rays show indeed very different yields and peak struc- 
tures. 

The (^He,a7) pick-up reaction reveals a cross section 
dominated by high £ neutron transfer. Here, the direct 
population of the residual nucleus takes place through 
one-particle-one-hole components of the wave function. 
Such configurations are not eigenstates of the nucleus, 
but they are rather distributed over virtually all eigen- 
states in the neighboring excitation-energy region. Thus, 
the neutron-hole strength for single particle levels away 
from the Fermi energy is distributed over a rather large 
range of background states. 

On the other hand the inelastic scattering (^IIe,'^IIe'7) 
reaction is known to populate mainly collective excita- 
tions with a slightly lower spin window. Collective exci- 
tations built on the ground state give rise to rather pure 
eigenfunctions and their strength is less spread over other 
eigenfunctions of the nucleus in the neighboring excita- 
tion energy region. 

In order to test if the number of 7-rays per cascade 
depends on the two types of reactions, we have evaluated 
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the average 7-ray multiplicity 



E 



(1) 



as a function of excitation energy E. The average 7-ray 
energy (E^) is calculated from 7 spectra selected at a 
certain energy E. 

Figure [3 shows the 7-ray multiplicity versus excitation 
energy. In spite of the different reaction mechanisms, 
the two reactions give similar results. In particular, the 
multiplicities (solid and dashed lines) of ^^Mo and ^^Mo 
are equal within their error bars, which gives support to 
the applicability of the Oslo method for both reactions. 

The experimental extraction procedure and assump- 
tions of the Oslo method are given in Refs. 0111 and 
references therein. The first generation (or primary) 7- 
ray matrix which is obtained as described above can be 
factorized according to the Brink- Axel hypothesis [30ll3l| 
as 



P{E, E^) cx p{E - E^)T{E^ 



(2) 



where p is the level density and T is the radiative trans- 
mission coefficient. 

The p and T functions can be determined by an itera- 
tive procedure 1^ through the adjustment of each data 
points of these two functions until a global minimum 
with the experimental P{E^ E-y) matrix is reached. It has 
been shown |0| that if one solution for the multiplica- 
tive functions p and T is known, one may construct an 
infinite number of other functions, which give identical 
fits to the P matrix by 

p{E-Ey) = Aexp[a{E - Ey)]p{E - Ey), (3) 
f{E-y) = Bexp{aE^)T{E^). (4) 

Consequently, neither the slope nor the absolute values 
of the two functions can be obtained through the fitting 
procedure. Thus, the parameters a, A and B remain to 
be determined. 

The parameters A and a can be determined by normal- 
izing the level density to the number of known discrete 
levels at low excitation energy and to the level den- 
sity estimated from neutron-resonance spacing data at 
the neutron binding energy E = i?„ |33j . The procedure 
for extracting the total level density p from the resonance 
energy spacing D is described in Ref. Q. Since our ex- 
perimental level-density data points reach up to an exci- 
tation energy of only E ~ i?„ — 1 MeV, we extrapolate 
with the back-shifted Fermi-gas model [SJ, |33 



Pbsfg{E) = T] 



12\/2ai/4[/5/4^/ 
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where a constant ij is introduced to adjust /Obsfg to the 
experimental level density at i?„. The intrinsic excita- 
tion energy is estimated by C/ — E — Ci — -Epair, where 
Ci = —6.6^4"" '^^ MeV and A are the back-shift param- 
eter and mass number, respectively. The pairing energy 



impair is based on pairing gap parameters Ap and A„ eval- 
uated from even-odd mass differenc es ISdl following the 
prescription of Dobaczewski et al. The level den- 

sity parameter is given by a = 0.21A°-®^ MeV~^. The 
spin-cutoff parameter tr/ is given by erf — O.OSSSaTA^/'^, 
where the nuclear temperature is given by 



(6) 



In cases where the intrinsic excitation energy U becomes 
negative, we set C/ = 0, T = and aj = 1. 

Figure 131 demonstrates the level density normalization 
procedure for the ^^Mo case. The experimental data 
points are normalized according to ^ by adjusting the 
parameters A and a such that a least fit is obtained in 
between the arrows. For the lower excitation region (see 
upper panel), one should take care only to include a fit 
region where it is likely that (almost) all levels are known. 
In practice, this means that the level density should not 
exceed ~ 50 levels per MeV. At the higher excitation re- 
gion (lower panel), the Fermi-gas extrapolation connects 
seamlessly the highest-energy data points with the level- 
density value determined from neutron-resonance spacing 
at Bn- Generally, the resulting normalization is not very 
dependent on the choice of the theoretical extrapolation 
function (Fermi gas) or the chosen fit region (~ 4.5 to 
- 7 MeV). 

Unfortunately, no information exists on the level den- 
sity at E = Bn for ^'^Mo. Therefore, we estimate this 
value from a systematics of other Mo isotopes where in- 
formation on the level density at i?„ exists. In Fig. 0| 
we plot all the known data points from the Mo isotopic 
chain. The odd- and even-mass molybdenum nuclei fall 
into two groups, both showing a decreasing level density 
as function of excitation energy. This behavior is rather 
counterintuitive since in a given nucleus the level density 
increases exponentially with excitation energy, and for 
neighboring nuclei one would naively expect quite simi- 
lar level-density curves. Two effects combine to result in 
the negative slope of the data points: (i) the decrease of 
single-particle level density when approaching the N = 
50 shell gap resulting in a decrease of the level density in 
general, and (ii) the increase of the neutron-separation 
energy with decreasing neutron number. For the nega- 
tive slope to emerge, both effects have to be rather pre- 
cisely of the same size for each step along the Mo isotopic 
chain. We have found no good physical explanation for 
this to happen, but we employ this fortuitous coincidence 
to estimate p{Bn) = (6.2±1.0)104 MeV'^ for ^^Mo from 
our phenomenological systematics^. The splitting of data 
points between even and odd Mo isotopes must not be in- 
terpreted solely as due to the effect of the pairing-energy 
shift of the level-density curves; the difference in the mag- 



^ This val ue also agrees within a factor of two with the systematics 
of Ref. SI 
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nitude of B„ between neutron-odd and -even isotopes 
also affects the magnitude of this spUtting. 

III. LEVEL DENSITY AND FINE 
STRUCTURES OF THE ENTROPY 

The present knowledge on level density is concentrated 
in mainly two regions; the low excitation region up to 
~ 2 MeV, studied in detail using spectroscopy and count- 
ing of known, discrete levels [221 and the region around 
the neutron separation energy studied by experiments on 
neutron resonances ^3 ■ Almost nothing is known of the 
level density in between the above mentioned regions, 
but it is possible to determine quite reliably two anchor 
points of the level density. 

Figure |5l shows the extracted anchor points (filled data 
points) for nine molybdenum isotopes together with the 
level density deduced from known discrete levels (solid 
lines) . The upper anchor point is simply determined from 
neutron resonance data. The lower anchor point, which is 
the average value of three data points, is determined such 
that a straight line on a logarithmic plot, going through 
the upper anchor point, provides an upper bound of the 
level density distribution of known levels. The algorithm 
is iterative and treats all nuclei similarly to ensure that 
the results are comparable. The straight line connect- 
ing the lower and upper anchor points, is defined by the 
constant temperature formula 

p{E) = Ce^/^ (7) 

with = (lnp2 — ln/0i)/(-E2 ~ Ei) and C = 

pi exp(— iJi/r). Details are given in |41j|. Provided that 
all the levels are measured at the excitation energy of the 
lower anchor point, we find from the plots of Fig. [Slthat 
the temperature- like parameter t drops from 1.05 MeV 
for the spherical ^^Mo to about 0.72 MeV for the weU 
deformed ^"^Mo nucleus. This figure also illustrates the 
excitation energy where one would expect the appear- 
ance of missing levels in spectroscopic work, typically if 
the density of levels exceeds 50 MeV^^. 

The level densities p{E) extracted from the eight re- 
actions are displayed in Fig. The data have been nor- 
malized as prescribed above, and the parameters used for 
93-98]y[Q p ijg^gj TablcHI We find that the 

normalization constant rj drops by one order of magni- 
tude when going from deformed to spherical nuclei. This 
means that the spherical ^'^Mo has about ten times lower 
level density than predicted by a global Fermi-gas model. 
As mentioned earlier, this effect is one of the reasons for 
the negative slope of the data points in Fig. 0] 

Our experimental data interpolate between the previ- 
ously known lower anchor point at ~ 2 MeV and about 
1 MeV below the upper anchor point at ^ 7 MeV. For 
the energy interval between ~ 6 MeV and ~ 7 MeV, we 
rely on models [3 Ull' In spite of the fact, that the 
final extrapolation of the level density up to the nuclear 
binding energy is model dependent, this only affects the 



average slope of the level density and does not affect the 
fine structure. This enables us to observe fine structures 
in the level density which are thought to reflect the break- 
ing of individual pairs. In an earlier work, we showed how 
a simple single-particle plus pairing model can qualita- 
tively explain the emergence of such fine structures |2lj |. 
Moreover, we have in the past investigated how pairing 
correlations are weakened in the presence of already un- 
paired nucleons, but also how these unpaired nucleons 
around the Fermi energy can increase the cost in energy 
to break up further nucleon pairs due to the blocking ef- 
fect of the Pauli principle jij . Our goal in the present 
work is to obtain experimental values for the critical tem- 
perature of the pair-breaking process. On the way, we 
will also investigate some other thermodynamical prop- 
erties, in particular the entropy, when going from spheri- 
cal to deformed nuclei. The generalization of the concept 
of temperature for a small system is not straightforward 
and has been discussed extensively in the literature (see, 
e.g., Ref. and references therein). Traditionally, tem- 
perature is introduced in slightly different ways in the 
microcanonical statistical ensemble (as a property of the 
system itself) and in the canonical statistical ensemble 
(as imposed by a heat bath). The temperature-energy 
relations for rare-earth nuclei (the caloric curves) derived 
within the two statistical ensembles display in general a 
different behavior since the nuclei under discussion are 
essentially discrete systems [T^ . 

In order to avoid the shortcomings imposed by the 
above mentioned statistical ensembles, a new approach 
for the caloric curves based on the two-dimensional prob- 
ability distribution P{E,T) has been proposed [12 lisj. 
This approach bypasses the well-known problem of spu- 
rious structures such as negative temperatures and heat 
capacities in the microcanonical ensemble. On the other 
hand more structures in the new caloric curve are evi- 
dent than in the canonical caloric curve. However, this 
new method is still not well settled and we will defer the 
discussion of caloric curves to another occasion. 

Within the microcanonical ensemble the experimen- 
tally obtained level density corresponds to the partition 
function which is simply the multiplicity f2 of accessible 
states. Thus, the entropy S of the system within the 
energy interval E and E + 6 is determined by: 

S{E) = kB In n{E) (8) 

where n{E) = p(E)/po and the Boltzmann constant is 
set to unity (fcs = 1) for simplicity^. In order to ful- 
fill the third law of thermodynamics; namely S 
when T — > 0, the normalization denominator is set to 
po = 1.5 MeV~^. Entropy as a function of energy can 



^ More precisely, multiplicity Q(E) is proportional to 
p(-B) (2(J(_E)) -I- 1), where {J(E)) is the average spin of 
levels with excitation energy E. However, in the present work, 
we neglect the weak excitation-energy dependence of {J{E)). 
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be defined and measured for small and mesoscopic sys- 
tems as well as for large systems. However, fluctuations 
in level spacings which are typical for small systems will 
make the entropy sensitive to exactly how the energy in- 
terval between E and E + 6E is chosen. Thus, Eq. (O 
is only useful if fl{E) is a sufRciently smooth function, 
i.e., for the case that its first derivative exists. Small 
statistical fluctuations in the entropy S may give rise to 
large contributions to the temperature T, which is de- 
fined within the microcanonical ensemble as 

Figure shows the entropy deduced within the micro- 
canonical ensemble for ^^'^"^Mo (upper panel) and ^^'^^Mo 
(lower panel) . The entropy curve plotted on a linear scale 
is essentially identical to the level-density curve on a log- 
arithmic scale. In general, the most efficient way to cre- 
ate additional states in atomic nuclei is to break J = 
nucleon Cooper pairs from the core. The resulting two 
nucleons may thereby be thermally excited rather inde- 
pendently to the available single particle levels around 
the Fermi surface. We therefore interpret, e.g., the step- 
like increases in entropy around 2-3 MeV excitation en- 
ergy in Fig. [3 as due to the breaking of the first nucleon 
Cooper pair. 

The entropies of odd-mass nuclei are higher than those 
of their even-even neighbors, even when their mass num- 
bers are lower. It is interesting to compare entropies 
between neighboring nuclei. The difference in entropy 

^S(^E^ 5'odd— mass 5'cvcn— even (-^0) 

is assumed to be a measure for the single particle entropy. 
The entropies of the almost spherical ^^Mo and ^^Mo 
(lower panel of Fig. O follow each other closely above 
E ^ 2.5 MeV. Here, the odd valence nucleon behaves 
almost as a passive spectator. For ^^'^^Mo, we find AS* > 
for £' > 2.5 MeV. The deformed case, (lower panel of 
Fig. I?!) exhibits an entropy difference of AS" > 1. This is 
less than the value of AS* ^ 2 found for rare-earth nuclei 

[mil. 

These observations can be explained qualitatively by 
the fact that the single particle entropy depends on the 
number of single particle orbitals that are available for 
excitations at a certain temperature. For ^^'^^Mo at low 
energies, the single neutron outside the closed shell can 
only occupy the two ^^^d 57/2 orbitals giving an en- 
tropy of In 2 ~ 0.7. For the deformed nucleus ^^'^^Mo, 
symmetry breaking results in a splitting of these two sin- 
gle particle orbitals into seven Nilsson orbitals, giving a 
total entropy of In 7 ^ 1.9, i.e., about one unit more than 
for the ^•^'^''Mo case. In the rare-earth region strong de- 
formation and intruder orbitals from other shells result in 
an even higher single particle level density, giving rise to 
an even larger single particle entropy. Although our sim- 
ple argument somewhat overestimates the observed single 
particle entropies, it provides a satisfactory explanation 



for the differences between the single particle entropies 
in the different cases. 

The entropy in atomic nuclei at low energies does not 
simply scale with the total number of nucleons. In the 
presence of pairing correlations, i.e., away from closed 
shells, the entropy scales rather with the number of un- 
paired nucleons at a certain excitation energy. When 
pairing correlations cannot form due to the large single 
particle level spacings around closed shells, an unpaired 
nucleon will behave almost as a passive spectator without 
contributing significantly to the entropy of the system. 

Around 5.5 MeV excitation energy, a bump is observed 
in the entropy curves for the lighter ^^'^^Mo nuclei. In 
light of the discussion above, it is unlikely that such a 
bump can be interpreted as the breaking of a nucleon 
Cooper pair. We propose that this bump is related to 
the sudden onset of neutron excitations across the N = 
50 shell gap. Due to the generally higher level density 
and the onset of deformation in the heavier Mo isotopes, 
the opening of the gg/2 shell might be less significant, 
leading to the effect being spread out in energy. However, 
such an effect might become visible in the lighter ^^'^^Mo 
nuclei. This interpretation is supported by the fact that 
the large transfer peak at 5.5 MeV excitation energy in 
the particle spectrum of the ^^Mo('^He,a7)^^Mo reaction 
at 45 MeV beam energy (see Fig. has been shown in 
an experiment at the Yale University Enge splitpole to 
be dominated by high £ transfer, most likely £ — Ah p6|. 

IV. PHASE TRANSITIONS 
A. Model 

In this Section we will utilize a recently developed ther- 
modynamic model [41I li^ which allows the investi- 
gation and classification of the pairing phase transition. 
The model is based on the canonical ensemble theory 
where equilibrium is obtained at a certain given tem- 
perature T. It can describe level densities for mid-shell 
nuclei in the mass regions A = 58, 106, 162, and 234. 

The basic idea of the model is the assumption of a 
reservoir of nucleon pairs. These nucleon pairs can be 
broken and the unpaired nucleons are thermally scattered 
into an infinite, equidistant, doubly degenerated single- 
particle level scheme. The nucleon pairs in the reservoir 
do not interact with each other and are thought to occupy 
an infinitely degenerated ground state. The nucleons in 
the single particle level scheme do not interact with each 
other either, but they have to obey the Pauli principle. 

The essential parameters of the model are: the num- 
ber of pairs n in the reservoir at zero temperature, the 
spacing of the single-particle level scheme e = 30 MeV /A, 
and the energy necessary to break a nucleon pair 2A = 
24 MeV/VA. Quenching of pairing correlations is intro- 
duced in this model by reducing the energy required to 
break a nucleon pair in the presence of unpaired nucle- 
ons. We assume that for every already broken nucleon 
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pair, the energy to break a further nucleon pair is reduced 
by a factor of r = 0.56 which is suggested by theoreti- 
cal calculations In order to simulate the effect of 
the TV = 50 shell closure on the ^ < 98 isotopes, we de- 
part from the global systematics for e and replace it by 
£/ = ea{A = 98) /a{A < 98) using the experimentally de- 
duced a values of Ref . j^^l ■ We use the same parameters 
for both protons and neutrons, keeping the proton pairs 
fixed to seven because there are 14 more protons outside 
the Z = 28 shell closure. 

The total partition function is written as a product of 
proton (Zt^), neutron (Z^), rotation (Zrot), and vibra- 
tion (.^vib) partition functions where the parameters for 
the collective excitations are the rigid moment of inertia 
A^ig ~ 34 MeV and the energy of one-phonon vi- 

brations ?iWvib = 1.5 MeV taken from spectroscopic data 
[3^. Thermodynamical quantities of interest can be de- 
duced from the Helmholtz free energy defined as 



F(r) = -rin(z,z,z„t^vib). 



(11) 



This equation connects statistical mechanics and thermo- 
dynamics. Quantities such as entropy, average excitation 
energy, and heat capacity can be calculated by 



S{T) 
(EiT)) 
Cv{T) 



F + ST 
fd{E)\ 



(12) 

(13) 
(14) 



respectively. 

In Fig. |S1 the Helmholtz free energy F, entropy S, av- 
erage excitation energy (£'), and heat capacity Cy are 
shown as functions of temperature for even-even, odd 
and odd-odd systems in the ^^Mo mass region. The 
free energy F and the average excitation energy (E) are 
rather structureless as functions of temperature. The 
odd-even effects are small: The even-even, odd and odd- 
odd systems have different excitation energies at the same 
temperature, where the even-even system requires the 
highest (E) to be heated to some given temperature T. 
Around Tc ^ 0.9-1.1 MeV the nuclei are excited to their 
respective nucleon separation energies. 

The entropy S and heat capacity Cy are more sen- 
sitive to thermal changes. The entropy difference AS' 
between systems with A and A ± 1 is a useful quantity. 
At moderate temperatures, it is approximately extensive 
(additive) and represents the single-particle entropy as- 
sociated with the valence particle (or hole) 
upper right panel, we find, e.g., that the nucleon carries 
a single-particle entropy of AS ~ 2.0 at T ~ 0.4 MeV. 

The shape of the heat-capacity curve is related to the 
level density. Traditionally, level-density curves have 
been described by the two component model as proposed 
by Gilbert and Cameron [s^. Within this model, the 
low energetic part is a constant-temperature level den- 
sity and the high energetic part is a Fermi-gas model. 



It has been shown in Ref. that the inclusion of a 
constant-temperature part in the level-density formula 
creates a heat-capacity curve as function of tempera- 
ture with a pronounced S shape similar to that shown 
in Fig. |S| With our model parameters, the maximum of 
the local increase in the Cy curve takes place at about 
T ~ 0.9 MeV. This temperature compares well with the 
temperature determined in the microcanonical ensemble 
from Eq. giving a temperature of T ~ 0.9 MeV for 
^^Mo (see also Fig.EJ. 



B. Comparison with experimental data 

Our model is described within the canonical ensemble 
while experimental data refer to the microcanonical en- 
semble. There are two ways to compare our model with 
experiments. Details are given in 141]. In this work we 
will make use of the saddle-point approximation |50l | 



pm) = 



exp(S') 



TV2^ 



(15) 



which gives satisfactory results for the nuclear level den- 
sity |4iLl4a|. 

FigureElshows the theoretical level densities calculated 
using Eq. H15|l . The agreement with the anchor points 
and the experimental level densities for ^^'^*Mo isotopes 
is satisfactory. Some of the model parameters could be 
adjusted more precisely, however, in this work we have 
chosen to use parameters taken from systematics. 

In order to investigate the behavior of the pairing cor- 
relations when approaching a major shell gap, we com- 
pare the canonical Cy curves which are based on the 
Laplace transforms of the experimental level densities. 
The curves are plotted on Fig. ^| for even 94,96,98]y[Q 
(upper panel) and odd 93,95,97]yjQ (lower panel) nuclei. 
The Cv curves resemble washed-out step structures and 
show an S shape as a function of temperature quite sim- 
ilar to the model calculation on the lower right panel 
of Fig. IHl Because of the averaging performed by the 
Laplace transformation discrete transitions between the 
different quasi-particle regimes, as observed within the 
micro-canonical ensemble, are hidden. Only the phase 
transition related to the quenching of the pair correla- 
tions as a whole can be seen. Details are given in ■ 

The canonical heat-capacity curves show local en- 
hancements around T ^ 0.5-1.0 MeV. Such enhance- 
ments were predicted in the calculations of Fig. |H1 and 
they are expected to be larger in the even-mass nuclei 
compared to the odd-mass neighbors. The experimen- 
tal heat capacities show this feature for the ^^'^^Mo pair, 
and up to a certain extend for the ^^'^^Mo pair, but it is 
not very obvious for the ^'^'^^Mo pair, where ^^Mo shows 
a more pronounced enhancement than expected. Ap- 
proaching the TV = 50 closed shell, the local enhance- 
ments become less and less pronounced. The general 
behavior of pairing correlations is that at shell closure 
there are almost no pairing correlations and, as particles 
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are added, the pairing correlations increase. Therefore 
the signature of a transition from a 'paired phase' to an 
'unpaired phase' when approaching a major shell gap be- 
comes less and less pronounced. We should note that very 
recently an alternative interpretation has been given 
These authors find that the S-shape can be accounted for 
as an effect of the particle-number conservation, and it 
occurs even when assuming a constant gap in the BCS 
theory. 

From the Cy curves, we have extracted the critical 
temperature for the quenching of pair correlations. The 
critical temperatures have been obtained by a fit of the 
canonical heat capacity of a constant-temperature level- 
density model to the data for the first 600 keV in tem- 
perature. The algorithm and its sensitivity are discussed 
in Ref. • The values obtained are plotted in Fig. ^2 
there is a tendency for the critical temperature to be 

Shghtly higher for odd 93,95,97]y[o ^j^^^j^ ^j. ^^^^ 92,94,96 jyj^ 

nuclei, similar to the local enhancement of the heat- 
capacity curve in the model calculation (see the lower 
right panel of Fig. ^ which is observed at higher tem- 
peratures for odd-mass Mo isotopes. The higher critical 
temperature for odd-mass nuclei is due to the Pauli block- 
ing effect of the unpaired quasiparticle which increases 
the distance to the Fermi surface for low-lying orbitals 
with coupled pairs and thus, increases the cost in energy 
to break pairs into more quasiparticles. Incidentally, the 
critical temperature for the quenching of pairing corre- 
lations coincides (by construction) quite well with the 
temperature-like parameter r of Fig. |5l 

A discontinuity of the heat capacity in a macroscopic 
system indicates a second-order phase transition accord- 
ing to the Ehrenfest classification; this is observed in the 
transition of a paired Fermion system such as a low- 
temperature superconductor or superfiuid '^He to their 
normal phases. Thus, the experimentally observed local 
enhancement of the heat capacity is interpreted as a fin- 
gerprint of a phase transition from a phase with strong 
pairing correlations to a phase where the pairing correla- 
tions are quenched T^. Shell-model Monte-Carlo calcu- 
lations 7] have shown that the pairing-phase transition 
is strongly correlated with the suppression of neutron 
pairs with increasing temperature. It has also been ob- 
served that the reduction of the neutron-pair content of 
the wavefunction is much stronger in the even-even than 
in the odd-mass isotopes, giving rise to the more pro- 
nounced S shape in the canonical heat-capacity curves 
in the even-even nuclei. The same odd-even difference in 
the heat capacity is also observed experimentally between 
leiDy and ^^^Dy, and ^^^Yb and ^^^yb [l^. 



C. Entropy as function of neutron number 

In order to study entropy as a function of neutron num- 
ber, we compare the micro-canonical entropy obtained by 
the saddle-point approximation of Eq. H15|l to our experi- 
mental data. In Fig.^jthe data are plotted as a function 



of the neutron number N (left panel) and as a function of 
the number of available neutrons in the reservoir (right 
panel). Although only qualitative agreement is achieved, 
some simple conclusions can be drawn. 

For the isotopes under investigation in this work, we 
see that the entropy at 1 MeV in both panels increases 
moderately as a function of the number of particles. The 
entropy at 7 MeV increases more rapidly and this is corre- 
lated to the evolution of the temperature-like parameter 
T, see Fig. O Both theoretically and experimentally, the 
odd systems show AS* — LOfc^ higher entropy than their 
neighboring even-even systems. 

The slopes at 7 MeV in the left panel of Fig. [T^l give 
dS/dN = O.bkg- Thus, going from ^^Mo to ^^Mo the 
level density drops when approaching the iV = 50 shell 
gap by a factor of ~ 0.03. This mechanism is also re- 
flected in the rj parameter of Eq. which drops from 
0.87 for ^^Mo to 0.08 for ^^Mo. 

As we already mentioned, the less pronounced S shape 
shows that the pairing correlations decrease when ap- 
proaching the = 50 shell gap. At the same time, the 
critical temperature for the quenching of pair correlations 
increases, which is the opposite of what one might expect. 
This effect can be explained by the increase in single- 
particle level spacing when approaching the A^ = 50 shell 
gap. We have already seen in the discussion in the pre- 
vious section, that this increase, together with the weak- 
ening pairing correlations, which fail to push the nuclear 
ground state sufficiently down in energy, lead to a de- 
crease in single particle entropy, see Figs. 171 and 1121 

Therefore, the increase in critical temperature for the 
quenching of pairing correlations when approaching the 
A^ = 50 shell gap is due to the competition between the 
weakening pairing correlations and the increasing single- 
particle level spacing. Just as the weakened pairing cor- 
relations in odd nuclei cannot compensate for the effect 
of Pauli blocking on Tc, they cannot compensate for the 
effect of an increase in single-particle level spacing on Tc 
when approaching a major shell gap. 



V. CONCLUSIONS 

Levels in ^^"^^Mo in the excitation-energy region up 
to the neutron-separation energy were populated using 
(■^IIe,Q;7) and ('^He,'^IIe'7) reactions. The level densi- 
ties of ^^~^^Mo were determined from their correspond- 
ing primary 7-ray spectra. Within the canonical ensem- 
ble, thermodynamical observables were deduced from the 
level density; they display features consistent with signa- 
tures of a phase transition from a strongly pair-correlated 
phase to a phase without strong pairing correlations. 
This conclusion is supported by recent theoretical cal- 
culations within shell-model Monte Carlo simulations by 
Alhassid et al. 0, H, E^l , where it is shown that the ex- 
pectation value of the pair operator decreases strongly 
around the critical temperature. However, we would 
like to point out, that other interpretations are not ex- 
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eluded. Different mcclianisms governing the thermody- 
namic properties of odd and even systems were studied. 
A simple, recently developed model for the investigation 
and classification of the pairing phase transition in hot 
nuclei has been employed and qualitative agreement with 
experimental data achieved. Using the saddle-point ap- 
proximation the experimental level densities of even-even 
and odd-even systems are reproduced. Estimates for the 
critical temperature of the pairing-phase transition yield 
Tc ~ 0.7-1.0 MeV. 
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TABLE I: Parameters used for the back-shifted Fermi-gas level density. 



Nucleus 


E 


a 


Ci 




D 


oiB„) 


'/ 




(MeV) 


(MeV~^) 


(MeV) 


(MeV) 


(eV) 


(lO^MeV"^) 




s^Mo 


2.080 


11.33 


-1.521 


8.642 


75 


9.99 


0.87 


''^Mo 


0.995 


11.23 


-1.526 


6.821 


1050 


3.10 


0.65 


^«Mo 


2.138 


11.13 


-1.531 


9.154 


105 


7.18 


0.46 


s^Mo 


1.047 


11.03 


-1.537 


7.367 


1320 


2.50 


0.34 


s^Mo 


2.027 


10.93 


-1.542 


9.678 




6.20" 


0.25 


^^Mo 


0.899 


10.83 


-1.547 


8.067 


2700 


1.27 


0.08 



"Estimated from systematics, see text. 
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FIG. 1: Charged ejectile spectra for ®^ ®*Mo in coincidence with 7-rays, labelled by the product nuclei. The arrows indicate 
the neutron separation energy B„. 
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FIG. 2: 7-ray multiplicity evaluated by Eq. (0 versus excitation energy. The individual spectra are labelled by the product 
nuclei. Solid and dashed lines represent (^He,Qf) and (^He,^He') reactions, respectively. 
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FIG. 3: Normalization procedure of the experimental level density (data points) of ^^Mo. The data points between the arrows 
in the upper panel arc normalized to known levels at low excitation energy (histograms) . In the lower panel they are normalized 
to the level density at the neutron-sepaxation energy (open triangle) using a Fermi-gas extrapolation (line). 



14 



C 



'Mo 



99 



Mo s7^o 



10^- 



95 



Mo 



93 



Mo 



94 s 



Mo 



o 



I .... I 



I .... I .... I .... I 



5 5.5 6 6.5 7 7.5 8 8,5 9 9.5 10 

Neutron -separation energy [MeVl 



FIG. 4: Level densities at the neutron separation energy. The unknown level density of ^'*Mo (open circle) is estimated from 
the slope of the data points of the odd-mass molybdenum isotopes. 
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FIG. 5: Level density of nine molybdenum isotopes. The histograms represent levels from spectroscopy [39|. A straight line is 
drawn from these levels to the level density at the neutron-separation energy which is determined by average neutron resonance 
spacings. The line represents the constant-temperature level-density formula, see text. 
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FIG. 6: Normalized level densities for ^*Mo. The open and filled circles are data from the ("^HejQ) and (^He,^He') reactions, 
respectively. 
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FIG. 7: Experimental entropy for ®^'^*Mo (upper panel) and ^^'®*Mo (lower panel) as function of excitation energy E. 
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FIG. 8: Model calculation for nuclei around Mo. The four panels show the free energy F, the entropy S, the thermal 
excitation energy {E), and the heat capacity Cv as a function of temperature T. The arrow at Tm ~ 0.9 MeV indicates the 
local maximum of Cv where the pair-breaking process takes place in the even-even system. The same parameter set is used 
for even-even (solid lines), odd (dashed lines), and odd-odd systems (dash-dotted lines). 
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FIG. 9: Calculated level density of ^*Mo (solid line) and ^'^Mo (dashed line) as function of average excitation energy (E). The 
big open circles and squares are experimental level-density anchor points from Ref. jSJ. The small filled and open circles are 
experimental data points measured with the (^He,a) and (^He,^He') reactions, respectively for the two isotopes. 



20 




FIG. 10: Observed heat capacity as functions of temperature in the canonical ensemble for the even ^''■^^'^^Mo (upper panel) 
and odd ^^'^^'^''Mo (lower panel) nuclei. 
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FIG. 11: Critical temperature for the quenching of pair correlations for Mo isotopes. 
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FIG. 12: Entropy extracted at 1 and 7 MeV excitation energy as function of neutron number N (left panel) and number of 
available neutrons in the model (right panel) for odd-even (open circles) and even-even (filled circles) molybdenum isotopes. 



